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We present a novel functional for spin density functional theory aiming at the description of non- 
collinear magnetic structures. The construction of the functional employs the spin-spiral-wave state 
of the uniform electron gas as reference system. We show that the functional depends on transverse 
gradients of the spin magnetization, i.e. , that in contrast to the widely used local spin density 
approximation the functional is sensitive to local changes of the direction of the spin magnetization. 
As a consequence the exchange-correlation magnetic field is not parallel to the spin magnetization 
and a local spin-torque is present in the ground state of the Kohn-Sham system. As a proof-of- 
principle we apply the functional to a Chromium mono-layer in the non-coUinear 120°-Neel state. 
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Density functional theory (DFT) |Ij is presently the 
most widely used approach to determine the electronic 
structure of atoms, molecules and solids. Shortly after 
the original formulation by Hohenberg and Kohn ,2. in 
terms of the electronic density n (r) the theory was ex- 
tended to include also the spin magnetization m{r) as 
fundamental variable [Sj . Spin density functional theory 
(SDFT) applies to Hamiltonians of the form 



■H = r + v + s + w, 



(1) 



where in addition to the kinetic energy T, the potential 
energy V and the electron-electron interaction energy VV 
(usually the Coulomb repulsion between two electrons) 
an additional energy contribution B due to an external 
magnetic field B(r) is considered, i.e. , 



B 



Mb d'^r B{r) ■ rh{r) 



(2) 



The operator m{r) — $^(r) cr<I>(r), representing the spin 
magnetization, is defined in terms of the spinor field 
(r) (/)j^(r)^ and the vector of Pauli matri- 
ces CT. As always in DFTs the success of the the- 
ory hinges on the availability of accurate and physically 
sound approximations to the exchange-correlation (xc) 
energy Exc[n,Tn\ - a functional of n{r) and m(r) in the 
case of SDFT. 

A straight-forward extension of the local density ap- 
proximation (LDA) J\ to SDFT yields 



E: 



LSDA 



[n,m]= Jd^r n{r) e'^f^ {n{r) , m{r)) 



(3) 



the so-called local spin density approximation (LSDA), 
with TO(r) being the magnitude of m{r) and e""'^ the xc 
energy of a spin-polarized uniform electron gas (UEG). 
For coUinear magnetism, i.e. , m(r) pointing in the same 



direction everywhere in space, a plethora of function- 
als was derived (cf. Ref. 5) improving over the collinear 
LSDA, however, much less is known about constructing 
functionals for non- collinear magnetism, where the direc- 
tion of m{r) is allowed to vary freely in space. In fact, 
most applications of non-collinear SDFT up-to-date are 
based on the idea, pioneered by Kiibler et al. [6^ , to apply 
collinear functionals to non-collinear systems by evaluat- 
ing the functional in a local reference frame with the local 
z-axis determined by the direction of m{r). The LSDA, 
defined in Eq. (|3]) , employs a local reference frame intrin- 
sically, which can be seen by evaluating the correspond- 
ing xc magnetic field 
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LSDA 



(r) 



dm /iBTn,(r) 



(4) 

By construction B^:^^^ is always aligned with m. The 
same is true for generalized gradient approximations 
(GGAs) employing the aforementioned rotation to a lo- 
cal reference frame. In recent years attempts were made 
to extend GGAs and meta-GGAs to non-collinear sys- 
tems without invoking a local reference frame in order 
to produce a Bxcir) which is non-collinear w.r.t. m(r) 
[7]. Since collinear functionals are usually formulated in 
terms of n^{r) and ^^(t) (as opposed to n(r) and m(r)) 
and gradients thereof, these approaches require a pre- 
scription mapping the gradient of m(r) (which is a 3x 3- 
matrix for non-collinear systems) to gradients of n^{r) 
and n^{r). Yet another approach was to consider the 
variations of the direction of m(r) perturbatively [SI [9]. 
In this letter we show that the very idea of the LSDA 
can be extended in a non-perturbative way to yield a new 
functional for SDFT depending on transverse gradients. 

In the LSDA the spin polarized UEG is chosen as refer- 
ence system to determine the local xc energy. Note that 
the LSDA does not employ the ground-state xc energy 
of the UEG, but instead the minimal xc energy of the 



2 



UEG under the constraint that its spin magnetization is 
mo- The new functional is based on the idea to consider 
a reference system with a non-collinear spin magnetiza- 
tion. In close analogy to the LSDA the local xc energy 
is determined from the UEG constrained to be in the so- 
called spin-spiral- wave (SSW) state [TU]. The SSW state 
of the UEG is characterized by a constant density no and 
a spin magnetization of the form 



mo(r) 



^scos(q-r) 
mo I ssin(<j'-r) 



(5) 



with s = sm{6) and 6 is the azimuthal angle between the 
rotating part (in the x-y plane) and the constant part 
(parallel to z-axis). The xc energy of the SSW UEG 
depends on four parameters: no, mo, s and q = \q\. As 
we will see below it is possible to define local s{r) and 
q{r) in terms of transverse gradients of m(r) which leads 
to the definition of the SSW functional 

Elf-^in, m] = Jd'r n(r) e^^^(n(r) , m(r) , s(r) , g(r)) , 

(6) 

where e^^^ is the minimal xc energy of the UEG un- 
der the constraint that it is in the SSW state speci- 
fied by n(r), m(r), s(r) and q{r). It is important to 
realize that the LSDA is included in this definition in 
the limits s — ?> or 9 — > 0, i.e. , e^f^{n, m,s,q ^ 0) = 



etT^{n,m,s^ 0, g) 



_ ^ssw 



elf^in,m,s^O,q = 0) = 
£™'^(n, m). This can be emphasized by rewriting the 
SSW functional as 

E^J'^in, m] = jd\ n(r) e^f{n{r) , m(r)) 

X (1 + S^c{n{r) , m(r) , s(r) , q{r))) , (7) 
where we introduced the spin gradient extension (SGE) 



Sxc{n,m,s,q) = 



SSW 



{n,m,s,q) -e™'f(n,m) 



runif 



(n, m) 



(8) 



Before we discuss the explicit form of the local s(r) 
and q{r) we briefly discuss global, i.e. , spatially indepen- 
dent, rotations of the internal (spin) space. These rota- 
tions correspond to transforming <t{r) — >-W$(r), where 
U is an element of SU(2) (a rotation of the internal [spin] 
degree of freedom). Note that spatial vectors, e.g. the 
paramagnetic current Jpir), are invariant under such 
internal rotations whereas spin vectors as m{r) trans- 
form as m(r) Rm{r), with R being the 3x3 rota- 
tion matrix corresponding to U. Since the kinetic en- 
ergy T and the interaction energy VV are invariant un- 
der global rotations of the internal space it follows that 
Sa;c[n, m] = Exc[n, Rm]. Considering infinitesimal spin 
rotations one obtains the so-called zero-torque theorem 



(9) 



0= d'rm{r) x B,,{r) 



which was first derived by Capelle et al. [TT] via the equa- 
tion of motion for the spin magnetization. It states that 
Bxc cannot exert a net torque on m. 

A simple rule to follow in order to ensure that explicit 
functionals for SDFT obey the zero-torque theorem is to 
write ExclnjTn] in terms of proper scalars, i.e. , spin in- 
dices have to be contracted with spin indices and spatial 
indices with spatial indices. This implies that the deter- 
mination of the local xc energy in terms of strictly local 
densities is exhausted by n(r) and m(r). Hence the lo- 
cal s{r) and g(r) have to be determined from properly 
contracted gradients of m(r). 

Let us first look at D{r) = |V®m(r) which corre- 
sponds to the total first order change of m{r). It can be 
split into a longitudinal contribution Di^{r) and a trans- 
verse contribution D^ir), i.e. , 

D{r) = ^—ADi^ir) + D^ir)) , 



m'^{r) 

D^ir) = \m{r)-{V(g>m{r)) |^ 
Drir) = \m{r) x (V(g)m(r)) p. 



(10) 

(11) 
(12) 



SSW UEG the two contributions are D^sw 



^ssw 
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rn^s q . 



where the meaning of longitudinal and transverse is 
defined by the local direction of m(r). For the 

and 

Both contributions are constant in 
space for the SSW UEG and hence play a similar role 
as the density no and the magnitude of the spin mag- 
netization mo, i.e. , they locally characterize the state. 
jjSSW vanishes because the spin magnetization in the 
SSW UEG only rotates. I?t does not vanish but it only 
determines the combination sq. 

Accordingly we look at the second order variation 
d{r) = |V^m(r) p. Again, it can be analyzed w.r.t. lon- 
gitudinal and transverse contributions 



dir) 



1 



m2(r) 
|m(r) 



(dL(r)+dT(r)), 



(V2m(r)) 



dLir) 

dT{r) = |m(r) x (V^m(r)) 

For our reference system this yields = 
dSSW ^mUl- s^ 



(13) 

(14) 
(15) 
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tjIqS q and 
The change of m(r) to first or- 
der is perpendicular to m(r), but to second order m(r) 
also changes in the direction of m{r) which explains why 
dl^"" does not vanish for the SSW UEG. In fact we 



see that provides the 



information as D^"^ , 



meaning, sq to some power. Adopting the convention 
that we obtain the characteristic parameters for the lo- 
cal xc energy using as few derivatives as possible s(r) 
and q{r) are determined through DT{r) and d^ir) by 




D^ir) +m^{r)dT{r) 



m' 



{r)DAr) 



(16) 
(17) 
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This completes the definition of the SSW functional 
Eq. (|6| or equivalently the SGE to the LSDA Eq. Q. 
By definition (c.f. Eq. (16l) the local s(r) is between 



[0, 1] in accordance with being the sine of an azimuthal 
angle. Furthermore we have the following hierarchy in 
the dependence of the SGE, Eq. ([t]), on the transverse 
gradients: i) If DT{r) = 0, the SGE correction is zero, ii) 
If DT{r) ^ and dT{r) — 0, the SGE correction is ob- 
tained from a planar SSW (s = 1). iii) If both transverse 
gradients are non-zero the SGE correction is obtained 
from a general SSW. 

We proceed by evaluating the xc magnetic field from 
the SSW functional, 



The zero-torque theorem, Eq. ([9]), is fulfilled by con- 
struction, however the new terms in Bxc{t) are non- 
coUinear w.r.t. m(r), i.e. , they provide a local torque 



m(r)xB,,(r) ^0. 



(22) 



This local torque accounts for the fact that in SDFT 
the KS system does not reproduce the exact many-body 
currents (spin and charge) but only the exact many-body 
densities n{r) and m(r) (review Ref. flT for a detailed 
discussion). Comparing the equation of motion for m(r) 
in the many-body system and the KS system one obtains 



V-J,,(r) = -2//Bm(r)xB,e(r): 



(23) 



Bg-+Bil, (18) 



where we split Bff^(r) into contributions coming from 
the dependence of on m, and dx, respectively. 

The explicit evaluation of Bff^(r) is straight-forward 
but rather lengthy. Here we will show the energetic con- 
tent in the Kohn-Sham (KS) system, i.e. , 



Ell 



-fiB Jd^r m{r)-Bxc{r) 
^dV n(r) {d„,elf^) mir) 
+ 4 Jd^r n{r) {dn^elt"^) D^{r) 
+ 4 /dV n{r) (a,,e^sw) ^^^^y 



(19) 
(20) 
(21) 



The first term (Eq. ( 19 )) is already present in the LSDA, 



whereas the other two terms (Eqs. (20),(21)) arise due to 
the inclusion of the SGE. 



16 



14 -\ 



12 



10 



— / 



, — --.-^x-v \ \ / -N \ \ I t ; / / 

-/ /. -W ) /^^ -N \ t / / 



I \ \ — 
I \ \ — 

■ 1 



6 -: 



"M / / -T / / 




/ /, — 




x\\ \ / — , 



^ I 
/ I 
^ I 
^ I 

1 • 

— t • 
-\ \ t 

\ \ t t 




- — / / 
- — / / 
- — / / 
-^--/ I 

till / / ^yy.'^ 
■ I t \ — -y^i^ 



1 1 / 
t 



I ^ 




\ 1 1 



-\ t 
— I 



\ X 



^ ^ ^ ^ 

\ i / / 

-xx\ I / -V 

W (-N \ 

' -/ ' ' 

-/ / 1 r 

. -/Ml 
w |-\ \ \ 
~xx\ I I -\ 

/ / /y'^-' 

- y 

\ \— ' 
l\ 

t I / / ^ ^ — 
\ \ I /z-^ — X 

X\\|/ -\\. 

XXXW \ \ \ t 

-/ML 




0.08 



0.07 



0.06 



0.05 



0.04 



0.03 



0.02 



0.01 



4 



10 



12 



14 



16 



FIG. 1. (Color online). Magnitude (heat map) and direction 
(arrows) of B^dr) for the Cr mono- layer (LSDA). 



where J^^^r) is the xc spin current defined as the dif- 
ference between the many-body spin current J(r) and 
the KS spin current Js(r). For the specific case of 
a vanishing external magnetic field (a case often stud- 
ied in SDFT), we know that the many-body currents 
vanish and hence the transverse part of B^c (r) deter- 
mines the irrotational part of the KS spin current |12) . 
In other words, a system exhibiting spontaneous non- 
coUinear magnetism is described by a KS system with 
non- vanishing irrotational spin current. 

The final step for a practical implementation of the 
SSW functional is the determination of the SGE from 
the SSW UEG. We obtained Sxc{n,m, s,q) using the 
random-phase approximation (RPA) for the SSW UEG. 
It is important to stress that we approximate the S^c 
with the RPA and not e^f^ . In this way the SSW func- 
tional reduces to the usual LSDA and not to the RPA- 
LSDA, i.e. , the LSDA parametrized using the RPA xc 
energy of the UEG. From ^ 65000 data points in the four- 
dimensional domain of S^c we constructed a polynomial 
fit for Sxc- 
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FIG. 2. (Color online). Same as FIG. [T] for the SSW func- 
tional. Note the intra-atomic non-coUinearity of Bxc{r). 
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FIG. 3. (Color online). Heat map of the z-component of 
m(r) xBxc{r) around a Chromium atom in the mono-layer. 
A non-vanishing KS spin current for the non-coUinear 120°- 
Neel state is implied (cf. Eq. (23l). Note the three-fold rota- 



tion symmetry inherited from the symmetry of the Chromium 
mono- layer. The zero-torque theorem (cf. (|9|) may be inferred 
from the pattern of negative (blue) and positive (red) local 
torques. 



As a proof-of-principle we implemented the SSW func- 
tional in the ELK code [13] in order to investigate the 
Chromium mono- layer in the 120°-Neel state. In FIGs.[l] 
and [2] we plot the magnitude and direction of B^c in 
order to illustrate the qualitative difference between the 
LSDA and the SSW functional. While the local spin 
magnetization m(r) are similar for the LSDA and the 
SSW functional, B^c obtained via the SGE exhibits much 
more structure compared to the LSDA B^c- As a re- 
sult (cf. FIG. [3]) the local torque does not vanish and a 
ground-state spin current is present in the KS system. 
Since the SSW functional is not restricted to small q{r) 
it accounts for the intra-atomic non-coUinearity. 

In conclusion we proposed a novel functional for SDFT 
depending on the first and second order transverse gra- 
dients of m(r). We emphasize that this functional is for- 
mulated in terms of a SGE to the LSDA, which vanishes 
in the case of a coUinear system. GGAs are also cor- 
rections to the LSDA, hence it is conceivable to employ 
the two corrections simultaneously. Since GGAs are con- 
structed having coUinear systems in mind, one may argue 
that the longitudinal gradients Di^{r), di,{r) should en- 
ter in the GGA part. We are confident that the SGE will 
improve the ab-initio description of materials exhibiting 
a non-collinear magnetic structure, which recently came 
in focus due to the observation of a Skyrmion lattice in 
the itinerant-electron magnet MnSi |T1]. 

While the corrections to the part of the Bxcir) parallel 
to m(r) will adjust the energetics, the perpendicular part 
oi Bxcif) describes the xc corrections to the spin current. 



which in turn is crucial for ab-initio spin dynamics. We 
hope that the functional presented in this letter will pave 
the road to a better description of domain wall motion 
and spin wave propagation from first principles in the 
framework of time-dependent SDFT. 

In both aforementioned scenarios it is important to 
have a numerically accessible functional which, given 
currently available computing facilities, implies the use 
of semi-local functionals. We demonstrated that non- 
coUinearity can be included by a generalization of the 
reference system employed in the LSDA and hence the 
numerical accessibility of the LSDA is retained in the 
SSW functional making it the ideal candidate for large 
scale quantum simulations. 
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